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The paper presents a dynamic solution method for dynamic minimum parametric networks flow. The solution method 
solves the problem for a special parametric dynamic network with linear lower bound functions of a single parameter. 
Instead directly work on the original network, the method implements a labelling algorithm in the parametric dynamic 
residual network and uses quickest paths from the source node to the sink node in the time-space network along 
which repeatedly decreases the dynamic flow for a sequence of parameter values, in their increasing order. In each 
iteration, the algorithm computes both the minimum flow for a certain subinterval of the parameter values, and the 
new breakpoint for the maximum parametric dynamic flow value function. 
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1 Introduction 


Dynamic flow problems where networks structure changes depending on a scalar parameter A are widely 
used (see for example Aronson ( 1989| l) to model different network-structured, decision-making prob¬ 
lems over time. These types of problems are arising in various real applications such as communication 
networks, air/road traffic control, and production systems. Moreover, in many applications of graph al¬ 
gorithms, including communication networks, graphics, assembly planning, and scheduling, graphs are 
subject to discrete changes, such as additions or deletions of arcs or nodes. In the last decade there has 
been a growing interest for such dynamically changing graphs, and a whole body of algorithms and data 
structures for dynamic graphs has been discovered: Fathabadi et al.| ( [2bl5| l, |He et al.| ( |2OT5l ), |Nassir et al. 
( |2014| l, IRashidi and Tsan^ ( |2015| l or |Zheng et al.|P015| . Further on, the next section presents some 
basic discrete-time dynamic networks terminology and notations and Section [^introduces the parametric 
minimum flow over time problem. In Section]^ the algorithm for solving the parametric minimum flow 
in discrete dynamic networks is presented and Section [^ilustrates how the algorithmon work on a given 
dynamic network. 
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2 Discrete-time dynamic network 

A discrete dynamic network G = {N, A, T) is a directed graph with N being a set of | | = n nodes 

i, A being a set of |A| = m arcs a and T being the finite time horizon discretized into the set H = 
{0,1, • • • ,T}. An arc a € A from node i to node j is denoted by (z, j). Parallel, as well as opposite 
arcs are not allowed in graph G. The following time-dependent functions are associated with each arc 
Q = ihj) G A: the upper bound (capacity) function u(i,j]S), u : Ax H ^ representing the 
maximum amount of flow that can enter the arc (z, j) at time 9, the lower bound function £(i,j',9), 
£: Ax H —>■ i.e. the minimum amount of flow that must enter the arc (z, j) at time 9, and the transit 
time function 9), h : A x H Time is measured in discrete steps, so that if one unit of flow 

leaves node z at time 9 over the arc a = (z, j) that one unit of flow arrives at node j at time 9 + h{i,j ; 9), 
where h{i,j',9) is the transit time of the arc (z,j). The time horizon T represents the time limit until 
which the flow can travel in the network. The network has two special nodes: a source node s and a sink 
node t. 

2.1 Time-space network 

For a given discrete-time dynamic network, the time-space network is a static network constructed by 
expanding the original network in the time dimension, by considering a separate copy of every node 
i G N at every discrete time step 9 G H. A node-time pair (NTP) (z, 9) refers to a particular node i G N 
at a particular time step 9 G H, i.e., (z, 9) G N x H. 

Definition 1 The time-space network of the original dynamic network G is defined as follows: 

(a) N'^ := {{i,9)\i G N, 9 G H}; 

(b) A^ := {ag = {(i,9), (j,9 -\- h(i,j]9)))\0 <9 < T - h{i,j]9),{i,j) G A}; 

(c) u^{ag) := u(a;9) for ag G A^; 

(d) £'^(ag) := £(a]9) for ag G A^. 

For every arc (z, j) G A with traversal time h{i,j',9), capacity u(i,j;9) and lower bound £(i,j;9), 
the time-space network G^ contains the arcs ((z, 9), (j, 9 -f ^*))), 0 = 0,1, • • •, T — 9) with 

capacities zz(z, j; 9) and lower bounds £{i,j]9). A (discrete-time) dynamic path P is defined as a sequence 
of distinct, consecutively linked NTPs. 

2.2 Parametric dynamic network 

Definition 2 A discrete-time dynamic network G = (N, A, T) for which the lower bounds £{i, j; 9) of 
some arcs (i,j) G A are functions of a real parameter A is referred to as a parametric dynamic network 
and is denoted by G = (N, A, u, £, h, T). 

For a parametric dynamic network G, the parametric lower bound function £: A x H x A\ ^ 
associates to each arc (z, j) G A and for each of the parameter values A in an interval [0, A] the real 
number £{i,j\(t; A), referred to as the parametric lower bound of arc (z, j): 

£(i,j-,9;X)=£o{i,j;9) + X-C(i,j;9), Ae[0,A], 9 G H, (1) 

where C : A x H —>■ ft is a real valued function, associating to each arc (z, j) G A and for each 
time-value 9 G H the real number referred to as the parametric part of the lower bound of 

the arc (z,j). The nonnegative value £o{i,j'A) is the lower bound of the arc (z,j) for A = 0, i.e.. 
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i{i,j\6]0) = £o{i,j;6). For the problem to be correctly formulated, the lower bound function of every 
arc {i,j) G A must respect the condition u(i,j\9) > £{i,j-,9;X) for the entire interval of the parameter 
values, i.e., V(i,j) G A,y9 G H andVA G [0,A]. It follows that 6*) > £o{i,j]9) and the parametric 
part of the lower bounds 9) must satisfy the constraints: C{i,j; 9) < [u{i,j; 9) — fo(*! J! (^)]/A, 

y{i,j) G A. 

3 Parametric flow over time 

The parametric dynamic flow value function v : N x H x [0, A] ^ iR associates to each of the nodes 
i G N, at each time moment 9 G H, a real number v{i; 9; A) referred to as the value of node i at time 9, 
for each of the parameter A values. 

Definitions A feasible parametric flow over time f{i,j\9]X) in a parametric dynamic network G = 
(TV, A, u, £, h, T) is a function / : A x iJ x [0, A] —> 3?“'' that satisfies the following constraints for all 
A e [0, A].- 


^ f{i,j;9;X)- ^ ^ /(j,T;iT;A) 


r v{i;9-,X), 

^ 0 , 


i = s,t 

1 Sjt 


,\/i G N,y9 G H; 


( 2 ) 

iii,j;9;X) < fii,j;9;X) <u{ij;9), V(i,j) S A, \/9 G H; (3) 

/(t,j;0;A)=O, V(T,j)eA, y9 G [T - h{i,j;9) +1, T]; (4) 

^v{s;9;X) =-^v{t;9]X) =v{X); (5) 

e^H e^H 

where f{i, j; 9; A) determines the rate of flow (per time unit) entering arc {i,j) 9>ne 9, for the parameter 

value X,y9 G {0,1, • • •, T} and VA e [0, A]. 

Definition 4 The parametric minimum flow over time (PmFT) problem is to compute all minimum flows 
over time for every possible value of X: 


minimise v(X) = v{s;9;X), for all As[0,A], 

6(e{0,l,-,T} 


(6) 


underflow constraints (H-Q- 

Definition 5 For the minimum flow over time problem, the parametric time-dependent residual network 
with respect to a given feasible parametric flow over time f is defined as Q{f) '■= (N, A(^f),T), with 
Mf) ■= ^^(/) U A-{f), where 


■= {ihj)\ihj) y <T - h{i,j-,9) with /(i, j; 6»; A) - j; 6»; A) > 0}; (7) 

A~{f) ■■= {ii,j)\ij,i) G A,39 <T - h{j,i-9) with u(j,i;9)-f{j,i-9;X)>0}. (8) 
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The direct arcs {i,j) G in G{f) have same transit times h{i,j; 9) with those in G while the 

reverse arcs (i, j) S A~ (/) have negative transit times h{i, j]9+ h{j,i\dy) = —h{j, i; 9), with (j, i) G A 
and 0 < 6* + h(j, z; 9) < T. 

Definition 6 For the parametric minimum flow over time (PmFT) problem, the parametric residual ca¬ 
pacities of arcs (i,j) in the parametric time-dependent residual network G{f) are defined as follows: 


= f{i,j\9]X) - l{i,j\9\\), {i,j)GA, 0 < 6< +/i(z, j; 6*) < T; (9) 

f{i,j;9-\-h{j,i;9)]X)=u{jfl;9)-f{j,i;9;X), {j,i) G A, 0 < 9-\-h{jfl;9) < T. (10) 

Definition 7 Given a parametric flow over time f(i, j; 9; A), the parametric residual capacity f{P; 9; X) 
of a dynamic path P is the inner envelope of the parametric residual capacity functions 0; X) for 

all arcs (i,j) composing the path and for all parameter X values: 

fiP;9;X) = min{f(z, j;0; A)|(z, j) S P}. (11) 

Definition 8 The transit time t{P\ 9) of a dynamic path P is defined by: 

t{P)= ^ h{i,j-9). (12) 


A dynamic path P is referred to as the quickest dynamic path if t(P) < t(P') for all dynamic paths 
P' in the parametric time-dependent residual network G{f ). 


4 Parametric minimum dynamic flow 

The approaches for solving the maximum parametric flow over time problem via applying classical al¬ 
gorithms can be grouped in two main categories: by applying a classical parametric flow algorithm (see 
Hamacher and Foulds (19891 and Ruhe ( 1991[)) in t he static time-space network (for this approach see the 
algorithm presented in Parpalea and Ciurea ( 2013|l) or by applying a non-parametric maximum dynamic 
flow algorithm (see the algorithms presented in Parpalea and Ciurea ( |2011| l) in dynamic residual networks 
generated by partitioning the interval of the parameter values (partitioning method presented by Avesalon 
jet al.| l. Given the powerful versatility of dynamic algorithms, it is not surprising that these algorithms and 
dynamic data structures are often more difficult to design and analyse than their static counterparts (see 
Orlin p013| l). Cai et al. (20071 proved that the complexity of finding a shortest dynamic flow augmenting 
path, by exploring the forward and reverse arcs successively, is 0{nmT^). For algorithms which explores 
the two sub-networks simultaneously, Miller-Hooks and Patterson] P004|l also reported a complexity of 


Nasrabadi and Hashemia 


(20101 


OlnfT^). By using special node addition and selection procedures, 
succeeded to reduce significantly the number of node time pair that needs to be visited. The worst-case 
complexity of their algorithm is 0(nT(n -\- T)). As far as we know, the problem of minimum flow over 
time in parametric networks has not been treated yet. 


4 .1 Parametric minimum dynamic fiow (PmDF) aigorithm 

The idea of the algorithm is that if the parametric residual capacities for all arcs in G{f) are maintained lin¬ 
ear functions of A, with no break points, the problem can be solved via a slightly modified non-parametric 






































Minimum parametric flow over time 


5 


algorithm. Firstly, if one exists, a feasible flow must be established. The most convenient is this to be 
done in the nonparametric network G' = {N, A, u, , h, T) obtained from the initial network by replacing 
the parametric lower bound functions with the non-parametric ones: = max{i{i, j', ; S', A)|A G 

[0,A]}, i.e. for < 0 and + A ' for 

L{i,j; 9) > 0. For finding a feasible flow 0) in G see the algorithms presented in Ahuja et al. 


( |1993| l, applied in the static time-space network. During the running of the algorithm, the subinter 
val of the parameter values is continuously narrowed so that the above restriction to remain valid. In 
each subinterval [Afc, Afc+i] of the parameter values, the linear parametric residual capacity of every arc 
{i,j) can be written as j',9\\) = ak{i,j',d) -f ■ (A — Afc) while the parametric flow is 

= fkihj',^) + Fk{i,j',S) • (A — Afe). As soon as the parametric time-dependent residual 
network G(/) contains no dynamic paths, the algorithm computes the minimum flow for the considered 
subinterval and then reiterates on the next subinterval of the parameter values, until A value is reached. 


(01) PmDF ALGORITHM; 

(02) BEGIN 

(03) find a feasible flow /(i, j ; 9) in network G; 
(04) SP:={0}; fc := 0; A*. := 0; 

(05) REPEAT 

(06) QDP(fc, Afc.SP); 

(07) A: := A: + 1; 

(08) UNTIL(Afc = A); 

(09) END. 


Algorithm 1: Parametric minimum dynamic flow (PmDF) algorithm. 


4.2 Quickest Dynamic Paths (QDP) procedure 

The successive shortest (quickest) dynamic paths procedure repeatedly performs the following operations: 

- Finds a shortest (quickest) dynamic path P in the time-dependent residual network G(/) based on the 
successor vector a \ N x H ^ (i,9) & N x H\ 

- Computes the parametric residual capacity r{P\ 6; A) of the dynamic path P; 

- Computes the first (in increasing order) value of the parameter A up to which the parametric residual 
capacity of the dynamic path remains linear without break points; 

- Updates the subinterval of the parameter values for which the minimum flow is determined; 

- Decreases the flow along the dynamic path and updates the time-dependent residual network. 

The algorithm ends when none of the source nodes, at any time moments 9 G {0,1, • • •, P}, is reachable 
from any of the sink nodes, i.e. there is no dynamic path from s to t. 


(01) PROCEDURE QDP(A:, Xu, BP)-, 

(02) BEGIN 

(03) FOR all 0 e IT DO 
(04) BEGIN 

(05) FORalH e AfDO(T(i,0) := (0,0); 

(06) FORall (i, j) e ADO 

0!k{i,j',9) ■= -Xk ■ C{i,j-,9)-, l3k(i,j-,9) := -C(i,j-,9)-, 

S + h(h 3 ',S)) ■■= - f{i,j-,9)-, 4- h{i,j-,e)) := 0 ; 

fk{i,j-,9) ■■= Fk{i,j-,e) -.= 0; 

(07) END; 

(08) C -.= 1; Afc+i = A; 

(09) LS(cr, G); 
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(10) WHILE(C = 1) DO 

(11) BEGIN _ 

(12) build P based on cr starting from (s, 6); 

(13) a := min{ofc(j, j;e)|(j, j) e P}\0 := m\n{j3k{i, j) S P, ak{i,j-,e) = a}; 

(14) {i,e):={s,e)- 

(15) WHILE(i ^ t) DO 

(16) BEGIN 

(17) (j,i7) := (t(*, 0); 

(18) < /3)THEN Afe+i := min{Afc+i, A*. + j;0) - o)/(/3 - /3fc(i, j; 0))}; 

(19) ak{i,j-,e) := ak{i,j;e) - a; /3fe(i,i;0) := lik{i,3\d) - /3; 

Ofe (i, *; 1?) := Ofc (i, *; 1^) + a; 

I3k{j,i\'&) ■= l3k{j,r,^) + /3; 

(20) iE((*,i) e Ai+(/))THEN/fe(i,i;0) ■.= fk{i, j-,e) - (x-Fk{i,r,e) ■.= Fk{i,r,e)-p-, 

(21) ELSE fk{j,i;^) ■■= fk{jp\’&) + a; Fk{jP\'&) ■= Fk{j,i\’^) + 

(22) (i,0) := 

(23) END; 

(24) EOR all 0 e -ff DO EOR alH e Af DO a{i, 0) ;= (0, 0); 

(25) LS((7, C); 

(26) END; 

(27) PP:=BPU{Afc+i}; 

(28) END. 


Algorithm 2: Quickest Dynamic Paths (QDP) procedure. 

4.3 Labels Setting (LS) procedure 

The label setting procedure uses transit time labels T{i, 0) associated to all nodes at each discrete time 
values. 

At any step of the algorithm, a label is permanent once it denotes the length of shortest augmenting 
path to a node-time pair, otherwise it is temporary. The LS procedure maintains a set L of candidate 
nodes in increasing order of their temporary labels, which initially includes only the sink nodes 
9 e {0,1, • • •, T}. For every node-time pair (j, tD) selected from the list, the arcs with positive residual 
capacity connecting (i, 9) to (j, d) are explored, where '9 = 9 + 9) if the arc connecting (i, 9) to 

(j, '9) is a forward arc and 0 < 0 = t? -f h{j, i; 9) if (i, j) is a reverse arc. At any iteration, the algorithm 
selects the node-time pair (z, 9) with the minimum temporary label, makes its transit time label permanent, 
checks optimality conditions and updates the labels accordingly. A transit time label r(z, 9) represents 
the length (transit time) of a shortest (quickest) dynamic path from (z, 9) if t(z, 9) < T{j, 9) + h{j, i; 9), 
y{i,j) € A{f). The process is repeated until there are no more candidate nodes in L. The transit time of 
the shortest (quickest) path P computed based on successor vector a is given by f. 

(01) PROCEDURE LS((7, C); 

(02) BEGIN 

(03) EORall0 e {0,1,---,T}DO 
(04) BEGIN 

(05) FOR all* e Af - I DO T{i, 0) := oo; 

(06) r(i, 0) := 0, L := L U {(i, 0)}; 

(07) END; 

(08) T := oo; 0 := oo 
(09) WHILE(L A 0) DO 

(10) BEGIN 

(11) select the first (j, 0) from L\ L := L — {{j, 0)}; 

(12) FOR alH e AT with {i,j) e A+(/) DO 



Minimum parametric flow over time 


1 


(13) FORall6»with6l + /i(i,j;6l) = i?DO 

(14) IF(r(j,i?) +/i(i, <T(i,e))THEN 

(15) BEGIN 

(16) T{i,e) := T{j,'d) + h{i,j-,0) 

(17) (T{i,e) := {j,0) 

(18) IF((i,6l) ^ L)THENL ■.= L U {{i, 9)}; 

(19) END; 

(20) FORalH e with (i, j) e A-(/) DO 

(21) IF)!? + h{j, v,'d) <T and T{j, 9) — h(j, i\ 9) < T(i, t? + h{j, i; t9))) THEN 

(22) BEGIN 

(23) T(i, 1? + h{j, i; 9)) := T{j, i9) — h{j ,«; 17); 

(24) i7 + h{j, i; i7)) := {j, i7); 

(25) IF((i,i? + /i(j,i;i7)) ^ L) THEN L := L U {(i, i7 + /i(j, i; i7))}; 

(26) END; 

(27) END; 

(28) f := mineg^o 1 ... r}{T(s,6»)}; 9 := mingg^o ^ ... r}{e|r(s, 0) = r}; 

(29) IF(t = oo) THEN C := 0; 

(30) END. 


Algorithm 3: Labels Setting (LS) procedure. 

Theorem 1 [Correctness of PmDF algorithm] PmDF algorithm computes correctly a parametric mini¬ 
mum flow over time for a given time horizon T and for X £ [0, A]. 


Proof: The partitioning type algorithm iterates on successive subintervals [Afc, Afc+i], starting with Aq = 0 
and ending with = A and consequently, the correctness of the algorithm obviously follows from 

the correctness of the Quickest Dynamic Paths (QDP) procedure which ends when none of the source 
node-time pairs is reachable from any of the sink node-time pairs, i.e. when there is no dynamic path from 
the source node to the sink node in the time-depending residual network. According to the classical flow 
decreasing path theorem presented by Ahuja et al. (19931, this means that the obtained flow is a minimum 
dynamic flow for the given time horizon. In fact, the algorithm ends with a set of linear parametric 
minimum flows and with the partition BP of the interval of the parameter values in their corresponding 
subintervals. □ 


Theorem 2 [Time complexity of PmDF algorithm] The parametric minimum dynamic flow (PmDF) al¬ 
gorithm runs in O(KnfmT^) time, where K -\- 1 is the number of X values in the set BP at the end of 
the algorithm. 

Proof: The building, as well as the updating of the time-dependent residual network requires an 0{mT) 
running time, since for each of the time values 0 < 0 < T all the m arcs must be examined. Labels Setting 
(LS) procedure investigates at most nT adjacent node-time pairs for each of the node-time pairs which 
are removed from the list L (i.e. 0{nT) times). Thus, the complexity of labels setting (LS) procedure 
is 0(n?T'^). Considering that in each of the iterations of the QDP procedure, for one of the time values, 
one arc is eliminated from the dynamic residual network, the algorithm end in at most 0{mT) iterations. 
On each of the iterations the procedure finds a quickest dynamic path with the complexity 0{nfT'^) and 
updates the time-dependent residual network in 0{mT) time. Thus, the total complexity of the quickest 
dynamic paths (QDP) procedure is 0{nfmT^). 

For each of the K subintervals [Afc, Afc+i], fc = 0,1, • • •, K—1 in which the interval [0, A] of the parameter 
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A values is partitioned in, the algorithm makes a call to procedure QDR Consequently, the complexity of 
the parametric minimum dynamic flow (PmDF) algorithm is 0{Kn^mT^). □ 


5 Example 

In the discrete-time dynamic network presented in Figure [T(a)| node 1 is the source node s and node 4 is 
the sink node f; the time horizon being set to T = 3, i.e. H = {0,1, 2,3}. 


10,00 


9,00 


8,00 


7,00 


6,00 


5,00 

0,00 0,20 0,40 0,60 0,80 1,00 

(b) The piecewise linear minimum flow over time 
value function for G. 

Fig. 1: Illustrating PmDF algorithm. 




(a) The discrete-time dynamic network G. 


For the interval of the parameter A values, set to [0,1], i.e. A = 1, the transit times h{i,j ; 6), parametric 
lower bound functions j; 0; A) = 9) + X ■ upper bounds (capacities) u{i,j; 6) and 

feasible flow / {i,j] 9) for all arcs in G are indicated in Tablefu 


Tab. 1: Characteristics of dynamic network G presented in Figure 1(a) 



eo{i,j-,e) C{i,j\S) u{ij;e) 

(1,2) 

1, 6» = 0 3, 0 = 0 -2, 0 = 0 5, 0 = 0 

2, 0>1 0, 0>1 0, 0>1 ^ 0, 0>1 

(1,3) 

1, O<0<2 1, O<0<2 5 2 0 -1 

i 0 ;; ' a 0 ;; 

(2,3) 

1, 0>O 0, 0>O 1 l = l 5 1’ l = l 

(2,4) 

1, O<0<2 2, 0 = 1 

2, 0 > 2 0, 0 > 0 0, 0 > 0 5 0 , 0 ^ 1 

(3,4) 

2, 0 < 0 < 2 J lZl<2 -2’ ^ 5 5 ’ 0 < 0 < 2 

0, 0^2 0<^_ 
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After the initialisation step, for k — 0 and for the corresponding initial value of the parameter Ao = 0, pro¬ 
cedure QDP is called for the first time. The successor vector is initialised for all nodes at all time values to 
cr{i, 0) := (0, 0) and the Labels setting (LS) procedure is called for finding a quickest dynamic path in the time- 
dependent residual network G{f). The distance labels are initialised to r(4, 6) 0, V0 G {0,1, 2, 3} and the 

set of candidate nodes is set to L := {(4, 0), (4,1), (4, 2), (4, 3)}. After setting transit time labels and finding 
successor values for all node-time pairs, the procedure computes the minimum label of the source node, r ;= 
min{r(l, 0), r(l, 1), r(l, 2), r(l, 3)} = min{2,1, oo, oo} = 1 with 0 := 1. Since r(4) 7 ^ 00 the procedure 
LS ends with the variable C keeping its initial value (7=1. Based on successor vector, the quickest dynamic path 
P := ((1,1), (3, 2), (2,1), (4, 2)) is built and its residual capacity r(P; 6^; A) = a -f (A — \k) ■ /? is computed with 
a = min{Q:o(l, 3; 1), ao(3, 2; 2), ao(2,4; 1)} = min{l, 2, 2} = 1 and /3 = —1. After the time-dependent residual 
network is updated and the parametric dynamic flow is decreased along the shortest dynamic path, the successor 
vector is reinitialised and procedure LS is called again. The next shortest dynamic path found by QDP procedure is 
P := ((1, 0), (2,1), (4, 2)) with r(P; A) = l-f A and both the time-dependent residual network and the parametric 
dynamic flow are accordingly updated. Then QDP procedure finds the new dynamic path P := ((1,0), (3,1), (4,3)) 
with a = 3 and j3 = Q. Since /lo(l,3;0) = —4 < /? = 0, the upper limit Xk+i of the subinterval of the pa¬ 
rameter is updated to Ai = min{l, (4 — 3)/(0 -f 4)} = 1/4. Finaly, after updating the time-dependent residual 
network and the parametric dynamic flow, based on successor vector found by LS procedure, QDP finds the last path 
P := (( 1 , 0 ), ( 2 , 1 ), (3, 2 ), (4, 3)) with a = 1 and /3 = 1 . The validity of the upper limit of the subinterval of the 
parameter values is tested and it is maintained unchanged since for the arc (2, 3) at 0 = 1, /lo(2, 3; 1) = —4 < /? = 1 
but Ai = min{l/4, (4 — 1)/(1 -|- 4)} = min{l/4, 3/5} = 1/4. At this point, after the updating step, no other dy¬ 
namic path can be found in the time-dependent residual network so that (7 := 0 is set, the breakpoints list is updated 
to BP = {0,1/4} and the first iteration ends with the minimum parametric flow over time /o(Li; d', A) computed 
for the subinterval [0,1/4] of the parameter A values. The PmDF algorithm increments variable k to the value fc = 1 
and a next iteration is performed. The evolution of the algorithm is presented in Table|^ The piecewise linear mini¬ 
mum flow over time value function for the discrete dynamic network (7, computed by Parametric minimum Dynamic 
Flow (PmDF) algorithm, is presented in Figure[T(b)| 


Tab. 2: The evolution of PmDF algorithm for the dynamic network (7 



k 

Afc 

P 

fk{P-,9-,X) 

A)c+i 




(1,1), (3, 2), (2,1), (4, 2) 

1 - A 

1 

Iteration 1 : 

0 

0 

(1,0), (2,1), (4, 2) 

1-f A 

1 




(1,0), (3,1), (4, 3) 

3 

1/4 




(1,0), (2,1), (3, 2), (4, 3) 

1-f A 

1/4 




(1,1), (3, 2), (2,1), (4, 2) 

1 - A 

1 

Iteration 2 : 

1 

1/4 

(1,0), (2,1), (4, 2) 

1-f A 

1 




(1,0), (3,1),(4,3) 

4-4A 

1 




(1,0), (2,1), (3, 2), (4, 3) 

l + \ 

3/5 




(1,1), (3, 2), (2,1), (4, 2) 

1 - A 

1 

Iteration 3 : 

2 

3/5 

(1,0), (2,1), (4, 2) 

1-f A 

1 




(1,0), (3,1), (4,3) 

4- 4A 

1 




(1,0), (2,1), (3, 2), (4, 3) 

4-4A 

1 
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